1 % 



BUNDfrSREPUB LIK D EUT^H 





REC'D 1 6 FEB 20M 



WIPO 



PCT 



Prioritatsbescheinigung uber die Einreichung 
einer Patentanmeldung 



Aktenzeichen: 
Anmeldetag: 
Anmelder/lnhaber: 
Bezeichnung: 



102 55 959.7 



29. November 2002 



Siemens Aktiengesellschaft, Munchen/DE 



Verfahren und Vorrichtung zur Analyse eines techni- 
schen Systems sowie Computerprogramm-Produkt 



IPC: 



G 05 B 17/00 




Die angehefteten Stucke sind eine richtige und genaue Wiedergabe der ur- 
sprunglichen Unterlagen dieser Patentanmeldung. 



MOnchen, den 26. November 2003 
Deutsches Patent- und Markenamt 
Der, President 



PRIORITY 
DOCUMENT 

SUBMITTED OR TRANSMITTED IN 
COMPLIANCE WITH RULE 17.1(a) OR (b) 



200219496 



Beschreibung 
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Verfahren und Vorrichtung zur Analyse eines technischen 
Systems sowie Computerprogramm-Erodukt 

Die Erfindung betrif ft ein Verfahren und eine Vorrichtung zum 
Entwurf eines technischen Systems sowie ein entsprechendes 
Computerprograrom-Produkt . 

Bei der Beschreibung eines technischen Systems, beispielswei- 
se einer Kraf twer ksanlage , werden verschiedene Parameter, wie 
z. B. Driicke, Massenflusse u.s.w. herangezogen. Die Parameter 
gehorchen bestimmten physikalischen Gesetzen, wie z. B. Mas- 
sen- Oder Energiebilanzen, welche sich durch ein Gleichungs- 
system ausdrticken lassen. Die L6sungen des Gleichungs systems 
sind die ZustandsgroBen des technischen Systems. Aus diesen 
Zustandsgr5Ben kSnnen wiederum fur den Betrieb des technx- 
schen Systems relevante DiagnosegroBen berechnet werden, wie 
z B. der Wirkungsgrad einer Kraf twer ksanlage. Der konkrete 
Zustand eines technischen Systems kann darOber hinaus durch 
Messungen erfasst werden. Die Messgr6Ben der Messungen kOnnen 
direkt den Wert einer ZustandsgrSBe wiedergeben; es konnen 
aber auch aus den ZustandsgrSBen abgeleitete MessgrSBen ge- 
messen werden. Beispielsweise kann die Temperatur eines tech- 
nischen Systems gemessen werden, wahrend die ZustandsgroBen 
des Systems die Enthalpie und der Druck sind. Zur Bestimmung 
der ZustandsgroBen aus den MessgroBen ftthrt man im allgemex- 
nen eine Messung durch und sucht die ZustandsgrSBen, welche 
das Gleichungssystem 16sen und deren abgeleiteten MessgrSBen 
meglichst nahe bei den durch die Messung ermittelten Messwer- 
ten liegen. Hierzu gibt es standardisierte Verfahren (siehe 
z.B. VDI-Richtlinie 2048). 

Es kann das Problem auftreten, dass aufgrund einer zu gerin- 
gen Anzahl von Gleichungen im Gleichungssystem bzw. einer zu 
geringen Anzahl von Messstellen einzelne ZustandsgrbBen oder 
einzelne DiagnosegrSBen unbestimmt bleiben. Ferner konnen 
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aufgrund von Messfehlern die ZustandsgrSfien bzw. die Diagno- 
segrSfien mit grofien Unsicherheiten behaftet. Es muss deshalb 
entschieden werden, durch welche Messungen die Genauigkeit 
bestimmter ZustandsgrSfien verbessert werden kann bzw. be- 
5 stimmte ZustandsgrSfien uberhaupt er£t bestimmt werden konnen. 
Hierzu wird ttblicherweise auf den Rat erfahrener Ingenieure 
zuruckgegrif fen, und die Vorschlage dieser Ingenieure lassen 
sich durch Simulationsprogramme iiberprufen. Hierzu sind je- 
doch zeitaufwendige Auswertungen erf order lich. 

10 - - 

Aufgabe der Erfindung ist es deshalb, ein Verfahren zum Ent 
wurf eines technischen Systems anzugeben, bei dem systema- 
tisch enrtittelt wird, wie die Messungen von einzelnen Mess- 
grSfien die Parameter des technischen Systems beeinf lussen. 

15 

Diese Aufgabe wird gemafi den Merkmalen der unabhangigen Pa- 
tentanspruche gelost. Weiterbildungen der Erfindung sind in 
den abhangigen Ansprtichen definiert. 

20 Das erfindungsgemafie Verfahren dient zum Entwurf eines tech- 
nischen Systems, das durch Parameter umfassend ZustandsgrSfien 
und von den ZustandsgrSfien abh&ngige DiagnosengrSfien charak- 
terisiert ist, Unter Entwurf wird hierbei insbesondere die 
Analyse und/oder Anderung des technischen Systems verstanden, 

25 insbesondere die Analyse und Anderung der im technischen Sysjg 
tern vorgenommenen Messungen. Das technische System wird dab«I 
durch ein Gleichun^ssystem beschrieben, wobei die Zustands- 
grSfien die LSsungen des Gleichungssystems sind. In den Ent- 
wurf des technischen Systems wird ein Messpark einbezogen, 

30 der erste MessgrSfien umfasst, wobei die ersten MessgrSfien im 
technischen System mit. einer vorgegebenen Genauigkeit gemes- 
sen werden, Ferner sind in dem technischen System zweite, von 
den ZustandsgrSfien abhangige MessgrSfien mit einer vorgegebe- 
nen Genauigkeit messbar. 

35 

In dem erf indungsgemSfien Verfahren werden erste Sensitivi- 
tatsgrSfien fur die ersten MessgrSfien und/oder zweite Sensiti- 
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vitatsgrCBen fur die zweiten Messgr6fien bestimmt, wobei zur 
Bestimmung der ersten SensitivitatsgrSBen ermittelt wird, in 
welchem Mafie eine Anderung der Genauigkeit der Messung der 
ersten MessgrSBen wenigstens einen ausgewahlten Parameter be- 
einflusst, und zur Bestimmung der zweiten SensitivitatsgrSBen 
ermittelt wird, in welchem MaBe die Messung der zweiten Mess- 
grSBen wenigstens einen ausgewahlten Parameter beeinflusst. 
Der Messpark wird daraufhin in Abhangigkeit von den ersten 
und/oder zweiten SensitivitatsgroBen derart verandert, dass 
die Genauigkeit einer oder mehrerer erster MessgrGBen veran- 
dert wird und/oder eine oder mehrere erste MessgrSBen aus dem 
Messpark herausgenommen werden und/oder eine oder mehrere 
zweite MessgroBen zum Messpark hinzugenommeri werden. Dieser 
veranderte Messpark wird zum Entwurf des technischen Systems 
eingesetzt . 

In einer bevorzugten Aus fuhrungs form wird vorzugsweise die 
Genauigkeit einer ersten MessgrSBe erhoht, wenn die erste 
SensitivitatsgrSBe dieser MessgrOfie in einem vorbestimmten 
Wertebereich liegt, und/oder eine erste MessgroBe wird aus 
dem Messpark herausgenommen, wenn die erste SensitivitatsgrS- 
Be dieser MessgrSBe in einem vorbestimmten Wertebereich 
liegt, und/oder eine zweite MessgrOBe wird zum Messpark hin- 
zugenommen, wenn die zweite SensitivitatsgrOBe dieser Mess- 
gr6Be in einem vorbestimmten Wertebereich liegt. Durch die 1 
Wahl von unterschiedlichen Wertebereichen kann somit das Ent- 
wurf sverfahren an unterschiedliche benutzerspezif ische Anfor- 
derungen in einfacher Art und Weise angepasst werden. 

In einer bevorzugten Aus fuhrungs form ist das technische Sys- 
tem durch ein Gleichungssystem H(x) = (HI (x) , Hn(x) = 0 

beschrieben, wobei x = (xl, ...xn) ein Vektor ist, der als 
Komponenten die ZustandsgroBen xi umfasst. An dieser Stelle 
sei angemerkt, dass alle im folgenden verwendeten Indizes i, 
j, k bzw. 1 nattirliche Zahlen darstellen. 
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Zur Durchftihrung des erf indungsgemafien Verfahrens werden in 
einer bevorzugten Ausf tihrungsf orm insbesondere folgande Mat- 
rizen berechnet: 

3 - eine Matrix N, die den Nullraum der Jakobimatrix von H auf- 

- eTnTtatrix W, so dass W '* die Inverse der Kovarianzmatrix 

der ersten Messgrbfien y^M 1st, wobei die Kovarianz- 
matrix alrf Eintrage die Kovarianzen cj= E((y 4 - E<y±>> <y 3 
0- - E (yj ))) aufweist,. wobei. E(y) der Erwartungswert von y 
ist • 

- eine Matrix M, welche die pseudoinverse Matrix zu A-lKDMr 
ist, wobei Db die Jakobimatrix der ersten Messgrbiien 
y±=b±(x) ist. 

Die Begriffe Nullraum, Jakobimatrix sowie inverse bzw pseu- 
doinverse Matrix sine, aus der Theorie der Matrixberechnungen 
bekannte Definitionen (siehe beispielsweise Gene H Golub, 
Charles F. van Loan: ..Matrix Computations", 3rd 
20 Baltimore, London; The Johns Hopkins university Press;. 1996) . 

in einer weiteren bevorzugten Ausfuhrungsform der Erfindung 
sind die in teohnisohen System bereohneten ersten Sensrtrvr- 

25 tatsgroAen jeweils das Verhaltnis der Genauigkeitsanderung 
eine! ausgewahlten Parameters zu der Genauigkeitsanderung e 
ner ersten MessgroBe, wobei der ausgewahlte Parameter erne 
ausgewahlte ZustandsgrOBe ist, die Ober die -sten Messgro.e 
bestimmbar ist. Das Verfahren zeiohnet sich hrerber daduroh 

30 aus, dass: 

- wenigstens einer der ausgewahlten Parameter eine ausge- 

* *. q _ .! q4- Hie tiber die ersten MessgroBen 
wahlte ZustandsgroBe ist, die uoei 

bestimmbar ist; 

35 - eine oder mehrere der ersten SensitivitatsgrSfcen O ^ De 
weils das Verhaltnis der Genauigkeitsanderung *m - 
AE(( Xl - E( X1 )) 2 )/xx der ausgewahlten ZustandsgrbBe Xl zu 
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der Genauigkeitsanderung A<Tjj 2 /yj = AE((yj - E(yj)) 2 )/yj ei- 
ner ersten MessgroBe yj darstellen; 

- die ersten SensitivitatsgroBen durch folgende Formel be- 
stimmt werden: 

V, -!j 

wobei r 2j das Element in der 1-ten Zeile und j-ten Spalte 
der Matrix N-M-W ist. 

In einer weiteren Ausfilhrungsf orm stellen die ersten Sensiti- 
vitatsgrSBen jeweils das Verhaltnis der Genauigkeitsanderung 
einer ausgewahlten DiagnosegrdBe zu der Genauigkeitsanderung 
einer ersten MessgrSBe dar, wobei die ausgewahlte Diagndse- 
grQBe iiber die ersten MessgroBen bestimmbar ist. Das Verfah- 
ren zeichnet sich hierbei dadurch aus, dass: 

- wenigstens einer der ausgewahlten Parameter eine ausge- 
wahlte DiagnosegrOBe ist, die ttber die ersten MessgrSBen 
bestimmbar ist; 

- eine Matrix Dd, welche die Jakobimatrix der DiagnosegrSBen 
dj.=di(x) ist, bestimmt wird; 

- eine Oder mehrere der ersten SensitivitatsgroBen <D y3 an 
jeweils das Verhaltnis der Genauigkeitsanderung Aonn 2 /d n = 
AE((d n - E(d n )) 2 )/d n der ausgewahlten DiagnosegrOBe d n zu 
der Genauigkeitsanderung Ao-jj 2 /yj = AE((yj - E(yj)) 2 )/yj ei- 
ner ersten MessgroBe yj darstellen; 

- die ersten SensitivitatsgroBen durch folgende Formel be- 
stimmt werden: 

wobei s Dj das Element in der n-ten Zeile und j-ten Spalte 
von bd-N'M-W ist. 

In einer weiteren bevorzugten Aus ftihrungs form stellen eine 
Oder mehrere der zweiten SensitivitatsgroBen jeweils die Va- 
rianz einer ausgewahlten ZustandsgrSfie bei der Hinzunahme ei- 
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ner zweiten MessgrSBe dar, wobei die ausgewShlte ZustandsgrS- 
fie uber die ersten MessgrSBen bestimmbar ist. Das. Verf ahren 
zeichnet sioh hierbei dadurch aus, dass: 

- wenigstens einer der ausgewahlten Parameter eine ausge- 
wahlte Zustandsgr5Be ist, die uber die ersten MessgroBen 
b£stimmbar ist; 

. eine Oder mehrere der zweiten SensitivitatsgrSBen D ewexls 
die Varianz <W der ausgewahlten ZustandsgraBe xx bei 
Hinzunahme einer zweiten MessgrOBe, deren Wert erne Zu- 
standsgroBe x, ist, mit der varianz <r, zu dem Messpark 

darstellen; , 
_ die zweiten SensitivitatsgrSBen durch folgende Formel be- 

stimmt werden: 

wobei m, die i-te Spalte der Matrix M *• N ist. 

in einer weiteren Ausftihrungsf orm der Erfindung stellen eine 
Oder mehrere der zweiten SensitivitatsgroBen ,ewexls die Va 
ri anz einer ausgewahlten DiagnosegrSBe bei Hinzunahme einer 
zweiten Messgr6Be dar, wobei die ausgewahlte 
Ober die ersten MessgrSBen beStimmbar ist. Das Verf ahren 
zeichnet sich hierbei dadurch aus, daSs: 

25 - wenigstens einer der ausgewahlten Parameter eine ausge- 
W 9 . MRo . ^ d , e ober die ersten MessgroBen 

wahlte DiagnosegroBe ist, die udsi 

bestimmbar ist; „,, or , 

. eine Matrix Dd, welche die Jakobimatrix der DiagnosegroBen 

d =d4 (x) ist, bestimmt wird; 
. eine Oder mehrere der zweiten SensitivitatsgroBen j swells 
die varianz o«->« 2 der ausgewahlten DiagnosegrSBe d. bei 
Hinzunahme einer zweiten MessgroBe, deren Wert erne Zu- 
standsgreBe x„ ist, mit der Varianz <z* zu dem Messpark 
35 darstellen; 
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- die zweiten SensitivitatsgrSBen durch folgende Formel be- 
st immt werden: 



wobei m± die i-te Spalte der Matrix M T - N T und q n die 
n-te Spalte der Matrix M r - N T - Dd T ist. 

Es kann nunmehr der Fall auftreten, dass der ausgewahlte Pa- 
rameter des technischen Systems eine ZustandsgrSfie ist, die 
nicht iiber die ersten Messgrofien bestiitimbar ist. In diesem 
Falle wird zunSchst eine zweite MessgrSBe ermittelt, deren 
Wert eine ZustandsgroBe ist und welche zu dem Messpark hinzu- 
zuftigen ist, damit der ausgewShlte Parameter eindeutig be- 
stimmbar ist. Das Verfahren, bei dem dieser Fall berucksich- 
tigt ist, zeichnet sich dadurch aus, dass: 

- wenigstens einer der ausgewahlten Parameter eine ausge- 
wShlte ZustandsgroBe ist, die nieht iiber die ersten Mess- 
grofien bestimmbar ist; 

- eine Matrix P, welche die orthogonale Projektion auf den 
Nullraum von A ist, bestimmt wird; 

- eine zweite MessgroBe ermittelt wird, deren Wert eine Zu- 
standsgrSfie x k ist und welche zu dem Messpark hinzuzuftigen 
ist, damit die ausgewahlte ZustandsgroBe eindeutig be- 
stimmbar ist; 

- eine der zweiten Sensitivitatsgrofien die Varianz a k -> x i 2 der 
ausgewShlten ZustandsgrSBe bei Hinzunahme der ermittelten 
zweiten MessgrSfie x k mit der Varianz er k zu dem Messpark 
darstellt; 

- die zweite Sensitivitatsgrofie durch folgende Formel be- 
stimmt wird: 




k 



l k 
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mit p = Pn lr wobei m die I-te Spalte der Matrix N * ist 
und m ± die i-te Spalte der Matrix M T - N T ist und p k die 
Jfc-te Spalte der Matrix P- N T ist. 

Es kann ferner der Fall auftreten, dass der ausgewahlte Para- 
meter eine DiagnosegrSBe ist, die nicht tiber die erste Mess- 
groBe bestimmbar ist. In diesem Fall wird zunachst eine zwei- 
te MessgroBe ermittelt, deren Wert eine ZustandsgroBe ist und 
' welche zu dem Messpark hinzuzufugen ist, damit die ausgewahl- 
• te DiagnosegrtSfie eindeutig best immbar ist.. Das Verfahren, . bex 
dem dieser Fall berucksichtigt ist, zeichnet sich dadurch 
aus, dass: 

- wenigstens einer der ausgewahlten Parameter eine ausge- 
wahlte DiagnosegrbBe ist, die nicht iiber die ersten Mess- 
groften bestiiranbar ist; 

- eine Matrix Dd, welche die Jakobimatrix der DiagnosegrSBen 

di=d±(x) ist, bestimmt wird; 

- eine Matrix P, welche die orthogonale Projektion auf den 
) Nullraum von A ist, bestimmt wird; 

- eine zweite MessgrSBen ermittelt wird, deren Wert eine Zu- 
standsgroBe x k ist und welche zu dem Messpark hinzuzufiigen 
ist, damit die ausgewahlte Diagnosegr6Be eindeutig be- 

stimmbar ist; 2 
25 - eine der zweien SensitivitatsgrSBen die Varianz a k - >dn der 
ausgewahlten Diagnosegr5Be d a bei Hinzunahme der ermittel 
ten zweiten MessgrSBe x k mit der Varianz <x* zu dem Mess- 

park darstellt; 
- die zweite SensitivitatsgrSfte durch folgende Formel be- 

30 stimmt wird: 

mit p - Pc„, wobei c n die n-ta Spalte der Matrix N *• Dd *, 
M k die Jr-te Spalte der Matrix M T - N T ist und p k die *-te 
35 Spalte der Matrix P* N T ist * 






9 

Bei den beiden zuletzt genannten ' Ausf tthrungsf ormen lauft die 
Ermittlung der zweiten MessgrOBe vorzugsweise derart ab, dass 
in der Matrix P ■ N T die Spalte gesucht wird, so dass p linear 
abhangig von dem durch diese Spalte dargestellten Vektor ist, 
5 wobei der Index der Spalte angibt, welcher zweite Messwert zu 
dem Messpark hinzuzuf tigen ist. 

In einer besonderes bevorzugten Ausgestaltung der Erfindung 
werden bei den Ausf uhrungsf ormen , bei denen als, zweite Sensi- 
10 tivitatsgroflen die Varianzen von ausgewahlten Parametern bei 
der Hinzunahme von zweiten Messgr6fien bestimmt werden f als 

•Standardabweichurig fur die zweiten hinzuzunehmenden Messgro- 
fien 1% des Wertes der zweiten MessgrGJJe angenommen. 

15 Es sei an dieser Stelle angemerkt f dass mathematisch die 

Richtigkeit aller im vorangegangenen verwendeten Formeln be- 
wiesen werden kann. 

Neben dem oben beschriebenen Verfaliren betrifft die Erfindung 
20 auch eine Vorrichtung zur Durchfuhrung des erf indungsgemSJJen 
Verfahrens. DarOber hinaus umfasst die Erfindung ein Compu- 
terprogrammprodukt, welches ein Speichermedium aufweist, auf 
welchem ein COmputerprogramm gespeichert ist, das auf einem 
Rechner ablaufbar ist und mit dem das erf indungsgemafie Ver- 

* ^^^^ 

mm fahren durchfiihrbar ist. 

Ausftihrungsbeispiele der Erfindung werden nachfolgend anhand 
der Zeichnungen dargestellt und erlautert. 

30 Es zeigt: 

Fig.l den schematischen Aufbau eines technischen Systems , 
das mittels des erf indungsgemaiien Verfahrens analy- 
siert wird, 

35 

Fig. 2 eine Prozessoreinheit zur Durchf iihrung des erf in- 
dungsgemaiien Verfahrens . 
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Das in Fig, 1 gezeigte technische System betrifft ein Heiz- 
systera eines Kraftwerks mit zwei hintereinander geschaltetfcn 
Heizflachen 1 und 2, wobei an diesen Heizflachen ein Gasstrom 
5 G und ein Wasserstrom W in jeweils entgegengesetzten Richtun- 
gen vorbeistrSmen. 

Das techpische System wird durch folgende Zustandsgrofien cha- 
rakterisiert : 

Massenfluss des Washers beim Eintritt in das 
Heizsystem; 

Massenfluss des Nassers beim Austritt aus dem 
Heizsystem; 

spezifische Enthalpie des Wassers beim Eintritt 
in das Heizsystem; 

spezifische Enthalpie des Wassers zwischen den 
beiden Heizflachen 1 und 2; 

spezifische Enthalpie beim Austritt des Wassers 
aus dem Heizsystem; 

Massenfluss des Gases im Heizsystem; 

spezifische Enthalpie des Gases beim Eintritt in 
das Heizsystem; 

spezifische Enthalpie des Gases- zwischen den bei- 
den Heizflachen 1 und 2; 



10 

m W,ein 



m W,aus 

15 



hw 9 ein 

25 

m G 

30 

^G 9 Mitte 



35 



spezifischen Enthalpie des Gases beim Austritt 
aus dem Heizsystem. 
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Die Zustandsgrofien sind Variablen eines Gleichungssystems 
H(x) = 0, welches folgende physikal'ischen Bilanzgleichungen 

umfasst : 

Massenbilanz des Wassers im Heizsystem: 

•v. 

m W 9 ein~ m W 9 aus =° 

Enthalpiebilanz an der ersten Heizflache: 
m G * { h G,ein ~ h G,Mitte )~ m W,aus *(%,aztf - h W 9 Mitte) = °' 
Enthalpiebilanz an der zweiten Heizflache: 
m G * (hG 9 Mtte ~ h G 9 aus)~ m W 9 em ' [ h W 9 Mitte " h W 9 em ) = 0 • 



Es wird folgend^r Soll-Arbeitspunkt des technischen Systems 
betrachtet, wobei die nachf olgenden Werte der ZustandsgrSfien 
eine Lbsung des obigen Gleichungssystems darstellen: 



m W,em 


m W,aus 


hjV,em 




hw,aus 


niQ 


^G,em 




^G,aus 


100 


100 


200 


300 


400 


50 


1000 


800 


600 



Neben den oben genannten Zustandsgr6fien ist das technische 
System ferner durch einen Diagnosewert charakterisiert , der 
im vorliegenden Fall den relativen Warmeubertrag des durch - 
str6menden Gases darstellt. Der warmetibertrag W kann durch 
folgende Formel beschrieben werden: 

jjr _ ^G 9 ein ^^G 9 aus 

Folgende erste Messgrofien wercien in dem technischen System 
mit einer Standardabweichung von jeweils 1% beziiglich des je- 
weiligen Soll-Wertes gemessen: 
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Enthalpief luss des Wassers beim Eintritt in das Heizsystem: 

5 Massenfluss des Wassers beim Eintritt in das Heizsystem: 
m W 9 ein ' 

Enthalpiefluss des Wassers beim Austritt aus dem Heizsystem: 

10 - . . .. _ 

m W 9 aus •v 



Massenfluss des Gases : 




15 m G ; 

Enthalpiefluss des Gases beim Eintritt in das Heizsystem: 
m G 'hG>em • 

20 

Mit der oben genannten Formel fur den relativen Warmetibertrag 
W ergibt sich mit den Sollwerten 0^ = 0,4. 

Aufgrund der verwendeten standardabweichungen von 1% fiihrt 
25 die Messung des relativen Warmeubertrags zu einem Messwert 
von 0,4 mit einer Standardabweichung von 0,0098. 

In einer ersten Variante des erf indungsgemafien Verf ahrens 
wird als SensitivitatsgrSfie eines ersten Messwerts jeweils 
30 das VerhSltnis der GenauigkeitsSnderung des Diagnosewertes W 
zu der Genauigkeit^anderung des ersten Messwertes bestimmt. 

Es ergeben sich folgende Werte : 



35 



SensitivitatsgrSfie des Enthalpief lusses des Wassers beim Ein- 
tritt in das Heizsystem: 0,167. 
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SensitivitatsgroBe des Massenf lusses des Wassers beim Ein- 
tritt in das Heizsystem: 0,0. 

SensitivitatsgroBe des Enthalpief lusses des Wassers beim Aus- 
5 tritt aus dem Heizsystem: 0,667. 

Sensitivitatsgr6fie des Massenf lusses des Gases : 0,0; 

SensitivitatsgrSBe des Enthalpief lusses des Gases beim Ein- 
10 tritt in das Heizsystem: 0,167. 

Man erkennt, dass die SensitivitatsgroBe des Enthalpief lusses 
des Wassers beim Austritt aus dem Heizsystem den groBten Wert 
aufweist. Das bedeutet, dass eine Anderung der Genauigkeit 
15 der Messung des Enthalpief lusses des Wassers beim Austritt 
aus dem Heizsystem den gr6Bten Einfluss auf die Genauigkeit 
der Messung des relativen Warmeubertrags hat. Folglich fiihrt 
eine Verbesserung der Messgenauigkeit des Enthalpief lusses 
des Wassers beim Austritt aus dem Heizsystem am wirkungs- 
20 vollsten zu einer Verbesserung der Genauigkeit des Diagnose- 
wertes. Demgegenuber haben die Messungen der Massenf lusse ei- 
nen Sensitivitatswert von 0 und deshalb keine Auswirkungen 
auf die Genauigkeit des Diagnosewertes W . 

In einer weiteren Ausf tihrungsf orm des erf indungsgemaBen Ver- 
fahrens werden als SensitivitatsgroBen die Varianzen des re- 
lativen Warmeubertrags unter der Annahme berechnet, dass im 
technischen System eine ZustandsgroBe zu den ersten Messgro- 
Ben als zweite Messgrofie mit einer Standardabweichung von 1 % 
30 beztiglich des Sollwertes hinzugenommen wird. Im Folgenden 

sind die Standardabweichungen (Wurzeln aus den Varianzen) bei 
der Hinzunahme der einzelnen Zustandsgrofien angegeben: 




35 



Hinzunahme von mjfr^i 0,0098; 
Hinzunahme von ntjy 9aus : 0,0098; 
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Hinzunahme von : 0,0095; 



10 



Hinzunahme von h Wt(ms : 0,0086; 
5 Hinzunahme von mo' 0,0098; 
Hinzunahme von h Gyein ' 0,0095; 
Hinzunahme von h G ,aus : 0,0062. 

Eine Hinzunahme von und ho, MIU »ird nicht betrachtet, 

da diese ZustandsgrSBen durch die MesagroBen nicht eindeutrg 
bestiaunt aind. Es iat eraichtiich, dasa die 
Mesaung von ^ die kleinste standardabwerchung fur den 
5 lativen Warmeubertrag ergibt. Folgiioh wird die Measung von 
CZ - ^ MesS * ar * dSr erStSn MeSSUngSn ^"^enon.aen. 

in einer weiteren Auafahrungaf orm der Erfindung warden nun- 
m ehr ZuatandagrSBen betrachtet, die nicht eindeutrg durch dre 
.0 ersten MeaagrOBen dea techniachen Syatema beatrnant srnd. Ea 
handelt aich hierbei um die zustandsgrOBen ™» • 

Ea wird nunmehr in eine* ersten Schritt ermittelt welche 

25 Rechnung gemafc Anspruch 11 durchgefuhrt . 

Es ergibt sich, dass eine Messung von %, M «e oder eine Mes- 
sung von W schon jeweils beide Zustandsgr5Aen -d 
festlegen. Bei der Annahme einer 
30 von 1 % fur die Messung von V,M«* bzw - ^,M« e er S lbt sicft ' 

- bei der Messung von h w , Mitte eine Standardabweichung fur 

von 3,0 und eine Standardabweichung fur % >M « e von 

17,34; 

35 
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- bei der Messung von %,Aff/te eine Standardabweichung fttr 
fyp,Mitte von 9 /°6 und eine Standardabweichung von h GMitte 
8,0. 

Hieraus ist ersichtlich, dass zur genauen Bestimmung von 
h^ 9 Mitte bevorzugt auch h^^ Ue gemessen wird, wohingegen zu ei- 
ner genauen Bestimmung von hQj^f Ue bevorzugt auch hoMtoe als 
Messung zum Messpark hinzugenommen wird. 

Das oben beschriebene Verfahren ermoglicht eine systematische 
und schnelle Suche nach neuen Mess-Stellen, durch welche die 
Genauigkeit von ausgewahlten Zustandsgrofien bzw. Diagnosegro- 
lien verbessert werden kann. Es muss somit nicht mehr auf die 
Erfahrung von Ingenieuren zuruckgegrif f en werden, urn zu ent- 
scheiden, welche Messgr6fien zu einem Messpark hinzugenommen 
werden sollen bzw. welche Messgenauigkeiten bevorzugt verbes- 
sert werden sollten. 

In Fig. 2 ist eine Prozessoreinheit PRZE zur Durchfiihrung des 
erfindungsgemaBen Verfahrens dargestellt. Die Prozessorein- 
heit PRZE umfasst einen Prozessor CPU, einen Speicher MEM und 
eine Input/Output-Schnittstelle IOS, die uber ein Interface 
IFC auf unterschiedliche Art und Weise genutzt wird: Ober ei- 
ne Grafikschnittstelle wird eine Ausgabe auf einem Monitor 
MON sichtbar und/oder auf einem Drucker PRT ausgegeben. Eine 
Eingabe erfolgt liber eine Maus MAS oder eine Tastatur TAST. 
Auch verfugt die Prozessoreinh.eit PRZE uber einen Datenbus 
BUS, der die Verbindung von einem Speicher MEM, dem Prozessor 
CPU und der Input/Output-Schnittstelle IOS gewahrleistet . 
Weiterhin sind an den Datenbus BUS zusatzliche Komponenten 
anschliefibar, z.B. zusatzlicher Speicher, Datenspeicher 
(Festplatte) oder Scanner. 
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Patentansprtiche 

1. Verfahren zxaiu Entwurf eines technischen Systems, das 

durch Parameter umfassend ^ustandsgrdfien und von den Zu- 
standsgr5Ben abhangige DiagnosegrSBen charakterisiert 
ist : 

- bei dem das technische System durch ein Gleichungssys- 
tem beschrieben wird, wobei die ZustandsgrSBen Losun- 
10 " gen des Gleichungs systems sind; _ . 

- bei dem ein Messpark umfassend erste MessgrSBen ana- 
lysiert wird, wobei die ersten MessgroBen mit vorge- 
gebenen Genauigkeiten im technischen System gemessen 
werden und von den ZustandsgrSBen abhangen; 

- bei dem zweite, von den ZustandsgrdBen abhSngige 
MessgroBen im technischen System mit einer vorbe- 
stimmten Genauigkeit messbar sind; 

- bei dem erste Sensitivitatsgrofien far die ersten 
MessgroBen und/oder zweite SensitivitatsgrSBen fur 
die zweiten MessgroBen bestimmt werdeh; 

- wobei zur Bestiiratiung der ersten SensitivitatsgrSBen 
ermittelt wird, in welchem MaBe eine Anderung der Ge 
nauigkeit der Messung der ersten MessgrSBen wenig- 
stens einen ausgewahlten Parameter beeinflusst, und 
zur Bestimmung der zweiten SensitivitatsgroBen ermit- 
telt wird, in welchem MaBe die Messung der zweiten 
MessgrSBen wenigstens einen ausgewahlten Parameter 
beeinflusst; 

- bei dem der Messpark in Abhangigkeit von den ersten 
und/oder zweiten SensitivitatsgrSBen derart verandert 
wird, dass die Genauigkeit einer oder mehrerer erster 
MessgroBen verandert wird und/oder eine oder mehrere 
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erste MessgroBen aus dem Messpark herausgenommen wer- 
den und/oder eine otier mehrere zweite Messgrolien zum 
Messpark hinzugenoramen werden; 

- bei dem der verSnderte Messpark zum Entwurf des tech- 
nischen Systems eingesetzt wird. 

2. Verfahren nach Anspruch 1, bei dem die Genauigkeit einer 
ersten MessgroBe erhoht wird, wenn die erste Sensitivi- 
tatsgrGBe dieser MessgroBe in einem vorbestimmten Werte- 
bereich liegt, und/oder eine erste MessgroBe aus dem 
Messpark herausgenommen wird, wenn die erste Sensitivi- 
tatsgroBe dieser MessgroBe in einem vorbestimmten Werte- 
bereich liegt, und/oder eine zweite MessgroBe zum Mess- 
park hinzugenommen wird, wenn die zweite Sensitivitats- 
grofie dieser Messgr6Be in einem vorbestimmten Wertebe- 
reich liegt. 

3. Verfahren nach Anspruch 1 oder 2, bei dem 

das technische System durch ein Gleichungs system 

H (x) = (H x (x) , . . . . ,H m (x) ) =0 beschrieben wird, wobei 

x= (xx, • • • r x n ) ein Vektor ist, der als Komponenten die Zu- 

standsgrofien x± umfasst. 

4. Verfahren nach Anspruch 3, bei dem die folgenden Matrizen 
berechnet werden: 

eine Matrix N, die den Nullraum der Jakobimatrix von H 
auf spannt, 

- eine Matrix W, so dass W T *W die Inverse der Kovarianz- 
matrix der ersten MessgroBen y±=b± (x) ist, wobei die 
Kovarianzmatrix als Eintrage die Kovarianzen a±j 2 = 
E((y± - E(yi))(yj - E(y-j))) aufweist, wobei E(y) der 
Erwartungswert von y ist; 

- eine Matrix M, welche die pseudoinverse Matrix zu 
A=W-Db-N ist, wobei Db die Jakobimatrix der ersten 
MessgroBen y±= s b±(x) ist- 
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5. Verfahren nach Anspruch 4, bei dem 

- wenigstens einer der ausgewahlten Parameter eine aus- 
gewahlte Zustandsgr6Be ist, die liber die ersten Mess- 
grofien bestimmbar ist; 

- eine oder mehrere der ersten Sensitivitatsgrofien O yjx i 
jeweils das Verhaitnis der Genauigkeitsanderung 
Ag u 2 /xi = AE((xi - E(xi)) 2 )/xi der ausgewahlten Zu- 
standsgrSfie xi zu der Genauigkeitsanderung Acrjj 2 /yj = 
AE ( (yj - E(yj) ) 2 ) /yj einer ersten MessgroBe yj darstel- 
len; 

- die ersten Sensitivitatsgrofien durch folgende Formel 
bestimmt werden: 

wobei rjy das Element in der I-ten Zeile und j-ten Spalte 
der Matrix N-M-W ist. 

6. Verfahren nach Anspruch 4 oder 5, bei dem 

- wenigstens einer der ausgewahlten Parameter eine aus- 
gewahlte Diagnosegrofie ist, die iiber die ersten Mess- 
grSfien bestimmbar ist; 

- eine Matrix Dd, welche die Jakobimatrix der Diagnose- 
groBen d±=di(x) ist, bestimmt wird; 

- eine oder mehrere der ersten Sensitivitatsgrofien O yj 
an jeweils das Verhaitnis der Genauigkeitsanderung 
AannVdn = AE ( (d n — E (d n ) ) 2 ) /d n der ausgewahlten Diagno 
segr6Be d n zu der Genauigkeitsanderung Acrjj 2 /yj - 
AE((yj - E(yj)) 2 )/yj einer ersten MessgroBe yj darstel 
len; 

- die ersten Sensitivitatsgrofien durch folgende Formel 
bestimmt werden: 
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wobei s n j das Element in der n-ten Zeile und j-ten Spalte 
von Dd-N-M-W 1st. 

7. Verfahren nach eineiq der Anspriiche 4 bis 6, bei dem 

- wenigstens einer der ausgewahlten Parameter eine aus- 
gewahlte Zustandsgrofie ist, die uber die ersten Mess- 
grofien bestimmbar ist; 

- eine Oder mehrere der zweiten SensitivitatsgrSfien je- 
weils die Varianz ajt-.> x i 2 der ausgewahlten Zustandsgrofie 
xi bei Hinzunahme einer zweiten Messgrofie, deren Wert 
eine Zustandsgrofie x k ist, mit der Varianz cr k zu dem 
Messpark darstellen; 

- die zweiten Sensitivitatsgrofien durch folgende Formel 
bestimmt werden: 



wobei m± die i-te Spalte der Matrix M r * N ist . 
8. Verfahren nach einem der Anspriiche 4 bis 7, bei dem 

- wenigstens einer der ausgewahlten Parameter eine aus- 
gewahlte Diagnosegrofie ist, die uber die ersten Mess- 
grofien bestimmbar ist; 

- eine Matrix Dd, welche die Jakobimatrix der Diagnose- 
grSfien d±=d±(x) ist, bestimmt wird; 

- eine oder mehrere der zweiten Sensitivitatsgr5fien je- 
weils die Varianz a k ->dn 2 der ausgewahlten DiagnosegroBe 
d n bei Hinzunahme einer zweiten Messgrofie, deren Wert 
eine Zustandsgrofie x k ist, mit der Varianz a k zu dem 
Messpark darstellen; 

- die zweiten Sensitivitatsgrofien durch folgende Formel 
bestimmt werden: 




200219496 



20 

wobei m± die i-te Spalte der Matrix M T - N T und q n die 
n-te Spalte der Matrix M r - N T > Dd T ist. 

5 9. Verfahren nach einem der Anspruche 4 bis 8, bei dem 

- wenigstens einer der ausgewahlten Parameter eine aus- 
gewahlte Zustandsgroiie ist, die nicht Uber die ersten 
MessgrSfien bestimmbar ist; 

10 - eine Matrix Py welche die orthogonale Projektion auf 

den Nullraum von A ist, bestimmt wird; 

- eine zweite Messgrofie ermittelt wird, deren Wert eine 
Zustandsgroiie x* ist und welche zu dem Messpark hinzu- 
zuftigen ist, damit die ausgew^hlte Zustandsgr&Jie ein- 

15 deutig bestimmbar ist; 

- eine der zweiten Sensitivitatsgr61ien die Varianz a k -> x i 2 
der ausgewahlten Zustandsgroiie bei Hinzunahme der er- 
raittelten zweiten MessgrcMie x k mit der Varianz <r k zu 
dem Messpark darstellt; 

20 - die zweite Sensitivitatsgrofie durch folgende Formel 

bestimmt wird: 

25 mit p = Pnx, wobei n 2 die 1-te Spalte der Matrix N T 

ist und m ± die i-te Spalte der Matrix M T - N T ist und 
p k die Jt-te Spalte der Matrix P- N T ist, 

10. Verfahren nach einem der Ansprijche 4 bis 9, bei dem 

30 

- wenigstens einer der ausgewahlten Parameter eine aus- 
gewahlte Diagnosegrofie ist, die nicht uber die ersten 
MessgrSBen bestimmbar ist; 

- eine Matrix Dd, welche die Jakobimatrix der Diagnose- 
35 groBen d ± =d± (x) ist, bestimmt wird; 
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- eine Matrix P, welche die orthogonale Projektion auf 
den Nullraum von A ist, bestimmt wird; 

- eine zweite MessgroBen ermittelt wird, deren Wert eine 
Zustandsgrdfie x k ist und welche zu dem Messpark hinzu- 

5 zufugen ist, damit die ausgewahlte Diagnosegrofie em- 

deutig bestimmbar ist; 2 
eine der zweiten SensitivitatsgrSBen die Varianz 0*-** 
der ausgewahlten Diagnosegrofie d n bei Hinzunahme der 
ermittelten zweiten MessgroBe x k mit der Varianz <r k zu 
10 dem Messpark darstellt; 

. die zweite SensitivitatsgrSBe durch folgende Formel 

• bestimmt wird: 
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mit p = PC, wobei c n die n-te Spalte der Matrix N *• 
Dd \ m k die Jc-te Spalte der Matrix M T - N x ist und p k 
die Jc-te Spalte der Matrix P- N r ist. 

11. Verfahren nach Anspruch 9 oder 10, bei dem in der Matrix 
P. W r die Spalte gesucht wird, sodass p von dieser Spalte 
linear abhangt, wobei der Index k dieser Spalte angibt, 
dass der zweite Messwert x k zum Messpark hinzuzunehmen 
ist, damit der ausgewahlte Parameter eindeutig bestimmbar 
ist. 

12. Verfahren nach einem der Anspruche 7 bis 11, bei dem 
die Standardabweichung <r* der zweiten MessgrSfie 1% des 
Wertes der zweiten MessgroBe ist. 

13. Vorrichtung zur Analyse eines technischen Systems, welche 
30 derart eingerichtet ist, dass ein Verfahren nach einem 

der vorhergehenden Anspruche durchfiihrbar ist. 

14. Computerprogrammprodukt, welches ein Speichermedium auf- 
weist, auf welchem ein Computerprogramm gespexchert ist, 

35 das auf einem Rechner ablaufbar ist und mit dem das Ver- 
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fahren nach einem der Ansprtiche 1 bis 12 durchfiihrbar 
1st. 
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Zusammenf assung 



Verfahren und Vorrichtung zur Analyse eines technischen Sys- 
tems sowie Computerprogramm-Produkt 

5 

Das erfindungsgemaBe Verfahren dient zum Entwurf eines tech- 
nischen Systems, das durch Parameter umfassend ZustandsgroBen 
und von den ZustandsgrOBen abhangige DiagnosengrOBen charak- 
terisiert ist. Das technische System wird dabei durch ein 
0 Gleichungssystem beschrieben, wobei die ZustandsgroBen die - 
LOsungen des Gleichungssystems sind. In den Entwurf des tech- 
nischen Systems wird ein Messpark einbezogen, der erste Mess- 
groBen umfasst, wobei die ersten MessgroBen im technischen 
System mit einer vorgegebenen Genauigkeit gemessen werden. 
L5 Ferner sind in dem technischen System zweite, von den Zu- 
standsgroBen abhangige MessgrOBen mit einer vorgegebenen Ge- 
nauigkeit messbar. In dem erf indungsgemaBen Verfahren werden 
erste SensitivitatsgroBen fur die ersten MessgroBen und/oder 
zweite SensitivitatsgroBen fur die zweiten MessgroBen be- 
20 stimmt, wobei zur Bestimmung der ersten SensitivitatsgroBen 

ermittelt wird, in welchem MaBe eine Anderung der Genauigkeit 
der Messung der ersten MessgrOBen wenigstens einen ausgewahl- 
ten Parameter beeinflusst, und zur Bestimmung der zweiten 
SensitivitatsgroBen ermittelt wird, in welchem MaBe die Mes- 
i% sung der zweiten MessgrOBen wenigstens einen ausgewahlten Pa- 
W rameter beeinflusst. Der Messpark wird daraufhin in Abhangig- 
keit von den ersten und/oder zweiten SensitivitatsgroBen der- 
art verandert, dass die Genauigkeit einer oder mehrerer ers- 
ter MessgroBen verandert wird und/oder eine oder mehrere ers- 
30 te MessgrOBen aus dem Messpark herausgenommen werden und/oder 
eine oder mehrere zweite MessgroBen zum Messpark hinzugenom- 
men werden. Dieser veranderte Messpark wird zum Entwurf des 
technischen Systems eingesetzt. 



35 Figur 1 
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